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Abstract. We study discrete Schrodinger operators {Hafiip){n) = ip{n — 1) -|- il){n -|- 1) -|- 
^ \ f{an + 9)'ip{n) on l'^{Z), where f(x) is a real analytic periodic function of period 1. We prove 

1-^ I a general theorem relating the measure of the spectrum of Ha^g to the measures of the spectra 

^ • of its canonical rational approximants under the condition that the Lyapunov exponents of 

Ha^g are positive. For the almost Mathieu operator (/(x) = 2Acos27ra:;) it follows that the 
measure of the spectrum is equal to 4|1 — |A|| for all real 0, X ^ ±1, and all irrational a. 

X 



1 Introduction 

Consider quasiperiodic operators acting on /^(Z) and given by: 

{H^,ei^){n)=ij{n-l)+ij{n + l) + f{an + 9)ij{n), n = ... ,-1,0,1, ... , (1.1) 

where f{x) is a real analytic periodic function of period 1. Denote by S{a, 9) the spectrum 
of Hafi. For rational a = p/q the spectrum consists of at most q intervals. Let S{a) = 
UeeM ^ ^'^^ ^) ■ ^ot^ ^hat for irrational a the spectrum of H (as a set) is independent of 6 
(see, e.g., [0]), and therefore S{a, 6) = S{a). In this paper we study continuity in a of S{a) 
and its measure. For sets, we will use | ■ | to denote the Lebesgue measure. 

The fact that various quantities could be easier to analyse and sometimes are even 
computable for periodic operators, Hp/^g, makes results on continuity in a particularly 
important. For example, the Aubry-Andre conjecture on the measure of the spectrum ||l[ 
states that for the almost Mathieu operator given by (|1 . 1|) with f{9) = 2Acos27r6', for 
irrational a and all real A, 6 there is an equality 

|5A(a,^)|=4|l-|A||. (1.2) 

Avron,van Mouche, Simon |p| proved that, for |A| 7^ 1, \Sx{Pn/<ln)\ -^ 4|1 — |A|| as g„ -^ cxd. 



and Last [|l3] established this fact for |A| = 1. Given these theorems, the proof of the 
Aubry-Andre conjecture was reduced to a continuity result. 

The continuity in a of S{a) was proven in [^ |^. Continuity of the measure of the 
spectrum is a more delicate issue, since, in particular, |5'(a)| can be discontinuous at rational 
a. Establishing continuity at irrational a requires quantitative estimates on the continuity 
of the spectrum. The first such result, namely the Holder- 1 continuity was proved in [0. 
That argument was improved to the Holder-1/2 continuity (for arbitrary / e C^) in ^ and 



subsequently used in [jT6|, O to establish (|1.2|) for the almost Mathieu operator for a with 



unbounded continuous fraction expansion, therefore proving the Aubry-Andre conjecture for 
a.e. (but not all) a. It was essentially argued in [^ that Holder continuity of any order larger 
than 1/2 would imply the desired continuity property of the measure of the spectrum for all 
a. Such continuity (more precisely, almost Lipschitz continuity, as in Theorem 3) was proved 
in ||13[ for the almost Mathieu operator with |A| > 14.5 (or dual regime), using exponential 
localization for Diophantine frequencies and delicate analysis of generalized eigenfunctions 
(positions of resonances and behavior between them). 

In the present paper we prove almost Lipschitz continuity of the spectrum for operators 
( |1.1|) with arbitrary analytic / under the assumption of positivity of the Lyapunov exponents. 



We would like to note that, unlike in |TB|], we do not rely on the theorem on exponential 



localization, and the present paper is essentially self-contained. In fact, the localization 
theorem for general analytic potentials, proved in [Q, establishes pure point spectrum for 
a.e. a without explicit Diophantine control, and therefore would not be applicable for our 
purposes. On the other hand, we would like to mention that, as the proof below shows, 
continuity of the spectrum requires simpler arguments than localization. 



Let Mn be the transfer-matrix of Hip = Eip : 

Then 

M„(9.i5) = (^-^(;" + '" -1). „ = .... -1.0.1..... 

Let Tn{9,E) = M.r,_i{e, E) ■ ■ ■ Mi{e, E)Mo{9, E) be the n-step transfer-matrix. The Lya- 
punov exponent '~f{E, a) of the family ( |1 . 1| ) is defined as follows: 

j(E,a)= lim - / In ||T„(^, E)||t/^ = inf - / \n\\TJ6,E)\\d9. (L3) 

n^oo n Jo n n Jo 

It is well defined by ( p.. 3D by Kingman's subadditive ergodic theorem (see, e.g., [^). Note 
that detM„(^,ii^) = 1 whence it follows that 'y{E,a) > 0. Our main result is 

Theorem 1. Let 'j{E, a) > for Lebesgue a.e. E. Then for any real 9, 






:l4) 



|^(a,^)| = lim 

71— >00 

where Pn/Qn is the sequence of canonical rational approximants to a. 

Remarks. 

1. Note that a.e. positivity of 7 implies that a is irrational. 

2. For a whose coefficients of the continued fraction expansion form an unbounded se- 
quence, the condition 'j{E, a) > is redundant. The result of Theorem 1 follows in this 
case (by a simple generalization of the argument in [|16|) from the Holder-1/2 continuity 
of the spectrum established in |Q. For a's whose continued fraction coefficients form a 
bounded sequence (henceforth, we denote the set of such a by fl) this continuity is not 
sufficient to imply Theorem 1. For such a we prove the theorem by first establishing 
the almost-Lipschitz continuity of the spectrum for which, in turn, we use positivity 
of the Lyapunov exponent. 

3. Note that it is not clear apriori that the limit above exists. 

4. While our proof can be easily adapted for arbitrary a we, for the reason above, con- 
centrate only on the case a G f2, which slightly simplifies certain arguments. 

We would like to add that recently a number of remarkable properties of quasiperiodic 
operators with analytic potentials have been established based on the positivity of the Lya- 
punov exponents [|T2| , ^, |^, ITof. Theorem 1, therefore, adds to the collection of general 
corollaries of positive Lyapunov exponents for analytic potentials. 



For the almost Mathieu operator Theorem 1 immediately implies 

Theorem 2. For every A, ^ G M, |A| 7^ 1, and every irrational a, the measure of the 
spectrum of the almost Mathieu operator H\^a,e is equal to 4|1 — |A||. 



This improves some results of [jT6|, |1^ and establishes the Aubry-Andre conjecture on 
the measure of the spectrum for all values of parameters in the noncritical case. For the 
critical case, A = 2, zero measure of the spectrum is known for the full measure set of a 



||TT| , 17|, however extending it to all irrational a remains an open problem. 

Proof: It was shown in that on a sequence of rational approximants Pn/Qn to any real 
a we have that \S\{pn/qn)\ (|A| 7^ 1) converges to 4|1 — |A||. Assume A > as the case A < 
is easily obtained by the transform 6-^6 + 1/2 and the case A = is trivial. We further 
assume A > 1 as the case < A < 1 can than be analyzed using duality (see, e.g. 0). For 
A > 1, it is known (e.g., ||^) that '~f{E, a) > for all real E and irrational a, which concludes 
the proof in view of Theorem 1. ■ 

In Section 2 we establish almost Lipschitz continuity of the spectrum and in Section 3 
we complete the proof of Theorem 1. 



2 Continuity of the spectrum. 

In what follows we shall often omit the indices a, 6, or E from the notation. Denote the 
basis in which H has the tridiagonal structure ( p..l| ) by {ei}'^_^. Let {E — i7)[a.^ .j.^] be the 
operator E — H restricted to the subspace spanned by {ci}^!^^. Set k = X2 + I — X1. Consider 
the polynomial in E: 

Pfc(«, e, E) = det((E - i/)[o,fc-i]), Po = 1. 

Let fk be the Fourier truncation of f{6) : 

h{e)= i: f,e^-'='. (2.5) 

j=-A;2 

As f{6) is an analytic function, its Fourier coefficients fj are of order e"'^-'. Hence |/(^) — 
fk{Q)\ < exp(— (ifc^), where d > depends on f{9), for k sufficiently large. Set z = exp{27ri9). 
Let Hk be the operator H with f{6) replaced by fk{0) and consider 

Pkia,z,E)=z'''det{{E-Hk)[o,k-i]). (2.6) 

Obviously, Pk{z) is a polynomial of degree 2k^ in z. Using induction for the three-term 
recurrence relation (which connects P„, P„_i, and Pn-2) and boundedness of f{9), we obtain 
for all z, \z\ = 1, and ii^ in a bounded set: 

\Pk\ - e-*' < \Pk\ < \Pk\ + e-*', \Pk\ < e^^ (2.7) 



where D < oo, dis somewhat decreased from the prior value but still positive, and k is larger 
than some ki. 

We shall now obtain the upper and lower bounds on \Pk{z,E)\ which we use later to 
establish the continuity properties of the spectrum. 

Lemma 1. (Estimate of \Pk{z,E)\ from below) For every aEQ, E&S, e'>0 
and at least one integer k out of each three consecutive integers s, s + 1, s + 2, s > ki, 
there exists C{e' , a) such that in any interval of length Ck^ there is an integer m such that 
|p^(g2^i(e+ma)^^)| > ^(^[E)-e')k^ rj.^^^ constant C IS independent of E. 

Remark. This statement, appropriately adjusted, holds under any Diophantine condi- 
tion on a. 

Proof: Since, as is easy to verify. 



we can rewrite ( |1.3|) in the form 



7(E) = inf^ Tin /i,,(^)d0, (21 



k k Jo 

where /ifc(^) = max{|Pfc(0)|, |Pfc-i(0)|, \Pk-ii9 + a)\, \Pk-2i9 + a)\}. Let Bk = {9 E [0, 1) : 
fik{9) > e(^(^)-^')'=}. Then, using (^ and {^, we obtain 

k^{E) < [' \nfik{9)d9 = f + f <\Bk\kD + {l- \Bk\){i{E) - e')k. 

Jo JBk J[0,l)\Bk 

Therefore, \Bk\ > e' /{D — ■y{E) + e') > 3ci{e') for all E in any bounded set (in particular, 
the spectrum of H.) Hence among each three consecutive indices s,s + l,s + 2 > ki, there 
is at least one, denote it k, such that 

\N, := {9 e [0, 1) : \Pk{e''''')\ > e^^)"^')'^}! > cie') > 

Denote by M the set of all such k. 

Since Pk is a polynomial of degree 2k^ in e^'^*^, the set Nk consists of less than 9k^ 
intervals. Let J be the interval in Nk with the maximum length. We have 

1^1 > ^. (2-9) 

It follows from the Diophantine properties of Q that 9 + ma (modi) G J for at least one 
m E I, where / is any interval of length Ck^. For reader's convenience we give a simple 
argument for the above. 

Let a G [0, 1) be an irrational and consider the continued fraction expansion a = 
[0; ai, 02, . . .], where a, are positive integers. Recall (e.g.,|T5[) that the convergents Pn/Qn = 
[0; ai, a2, . . . , a„] have the properties 

7)0 — 

a = — ^ , |(5| < 1, qn = anqn~i + qn-2, n = 3,4, .... (2.10) 

qn q-nqn+i 



and the greatest common divisor gcd{pn, q-n) = 1- 

For a fixed n, consider the set 9j = 9 + aj(mod 1), j = 0, 1, . . . , g„ — 1. Let 9j and 9j' be 



nearest neighbors. We have, using (|2.1CI|) , 



aj - 



JPn 



JS 



QnQn+l 



< 



J 



< 



1 



QriQu+l Qn 



Since the points ■^^(mod 1), j = 0, 1, . . . , g„ — 1 are the same as -^-(mod 1), j = 0, 1 



up to the ordering, we have, using ( 2.11| ): 



max \9j — 9j'\ < 



Qn 



J = 1,2, 



Qn 



.Qn 



(2.11) 



, ■ ■ ■ 1 Qn -L; 



(2.12) 



Therefore, if ci/9k^ > 3/g„, there exists < m < g„ — 1 such that 9m G J- Accordingly, 

choose n so that 

27k^ 
Qn-i < < Qn- (2.13) 

Cl 



For a G f2 we have that a„ < B(a) for all n. Then 

Qn = anQn-i + '?n~2 < 25(a)g„_i < Ck^. 



(2.14) 



where C = 60B{a)/ci. Therefore, at one of any g„ (as defined by (|2.13|) ) of the points 9m 
with consecutive indices we have the needed estimate for Pfc(e^'^*^). Inequality ( |2.14|) implies 
the statement of the Lemma. ■ 



Lemma 2. (Estimate of \Pk{z,E)\ from above) For any 6,6 > 0, there exists a set 
F{6,e) C S such that \F{6,e)\ < 6 and for any a E Q, E E S\F{6,e), and sufficiently large 
k (k > k2{5,a,e)) we have \Pk{z,E)\ < e^'^^^^^^^'' . The constant k2 is independent of z and 
E. 

Remark. Lemma 2 actually holds uniformly, for all ii^ in a bounded set, 0. However, 
this result requires a relatively complicated argument, while a nonuniform statement above 
follows immediately from a general theorem on uniquely ergodic dynamical systems and is 
sufficient for our purposes. 

Proof: It is proved in that for any continuous subadditive cocycle /„ on a uniquely 
ergodic system, limsup-/„(x) < lim - / /„(i/i(x) uniformly in x. Therefore, for large k we 
have \Pk{,z)\ < e(Tfc(^'°)+^'^'= for all |z| = 1. 

Let Dn = {E e S ■.'ik > n'i\z\ =1, \Pk{z,E)\ < 6^^)+")^}. Then we have 1 5 \ D„ | ^0 
as n ^ oo, which completes the proof of Lemma 2. ■ 



We are now ready to formulate the result about continuity of the spectrum. 



Theorem 3. Suppose a & Vt and '~f{E,a) > for Lebesgue a.e. E. Then for any 
6 > 0, there exists a set A{5) C S{a) such that \A{5)\ < 5 and for any E G S{a) \ A{5) and 
sufficiently small difference \a — a'\ < h{a, 6) we can find E' G S{a') such that 

\E' -E\ < c{a,5)\{a - a')\n^ \a - a'\\. (2.15) 

The constant c is independent of E, E' , a' . 

Remarks. 



1. The exponent 3 in ( p. 15 ) can be replaced by 2 + e with any e > 0, as the reader 



can easily verify starting with equation (|2.5|) . For the almost Mathieu case, 3 can be 
replaced by 1. 

2. A similar statement (with 3 replaced by a higher power) holds under a standard Dio- 
phantine condition on a. 

3. This statement can be made uniform in E' G S* if we require 7(-E) to be bounded away 
from (as in the case of the almost Mathieu operator) and use the uniform version of 
Lemma 2. 

Proof: As in P, ||, [l^, [T^, for a given E G 5(0;) we construct an approximate eigenvector 
for Ha- In order to obtain almost-Lipschitz continuity we need the error of approximation 
to be exponentially small, as in [|l^. 

Fix aeVt. Let G = {E G 5 : 7(E) < g} for some 5- > 0. Then IG] ^ as 5- ^ 0. Let g 
be such that |G| < 5/2. Choose positive e\ e" so that e' ^e" < g/lQ. Set A{6) = F(|, e") UG. 
Fix E e S\ A{S). Set K = {k e M,k> 2k2{5, a, e") + 3} with ^2 from Lemma 2, and M 
defined in the proof of Lemma 1. 

For any Eq E S let G{x,y) be the matrix elements of G = {{H — Eo)[xi,x2])^^- Using 
Cramer's rule and ( p7|) at the last step, we obtain: 



iG(x,xoi< '^--;L., 'T'!-1 , iG(x,x2)i< '"^^--^^ ''^°^'+' 



|Pfe(e2-^-i,Eo)|-e-*' ' i-v-.-^yi - |P;,(e2-^-i,Eo)|-e-'^'=' ' 

(2.16) 
where k = X2 — xi + 1 > ki. 

We shall now fix xi, k, and Eq. First take a k E K. By Lemma 2, |P[(,fc_|_j)/2](-2, -E)| < 
exp {'y{E) + e"){k/2 + 1) for j = —3, —2, . . . , 2, all z. We can find a neighborhood of E, 
rk{E) of diameter smaller than exp{—kg/4), so that for any E" G rk{E) we have for all z 
{\z\=iy. |Pp+,)/2](S")| < |P[(fc+,)/2](^)|exp(£"(fc/2 + l))and|Pfe(E")| > \PkiE)\exp{-e'k). 
Now take a generalized eigenvalue Eq G rk{E). Let ip be the corresponding generalized 
eigenfunction, that is a formal solution to the equation Hip = Eoip with |'0(x)| < C{\x\ + 1) 



for all X. By [p!8| , |^ the set of Eo which admit such solutions is dense in the spectrum. There 
holds: 

maxM4 = ^l:^ = ,^<oo (2.17) 

^ \x\ + l IXmaxI + 1 

for some point Xmax- We normalize ip so that R = 1. 

As is easily seen by considering {H — _Eo)[xi,x2]'^) we have for x G [xi,X2]: 

- iIj(x) = G{x, xi)ip{xi - 1) + G{x, X2)V'(x2 + 1). (2.18) 

Define the interval / = [xmax — k — Ck^, x^ax " k], where C = C{6', a) is the constant from 
Lemma 1. Then |/| = Ck^, and in view of Lemma 1, we can find Xi = m in I (which 
fixes position of the interval [xi, 0:2]) such that |Pfc(e^'^*^'", E)\ > exp (7(-E) — e')k and hence 
\Pk[e'^'^^^"^ , Eq)\ > exp {'y{E) — 2e')k. We now evaluate 'il){x) at the midpoint of the interval 
[a;i,a;2]: xq = xi + [{k — l)/2] and at its nearest neighbors. Using Lemma 1 to evaluate 
the denominator in ( p.l6| ), we obtain |G(xo,a;i)| < 2 exp ((— 7(_E) + e)k/2 + 7(-E) + e), e = 
4(5' + e") and the same estimate for G{xo^X2) for sufficiently large k E K (depending on d). 
Applying (|2.18|) , we get 

\i,{x, - 1)1, |^(a;o)| < 4e-(^-^)^/2+^+^(|a;^,J + (C + 2)k^) (2.19) 

for sufficiently large k E K. 

Now let /' = [xmax + k, Xmax + k + Ck^]. Similarly, we choose an interval [x-^^, Xg] so that 
at Xi G /' we have a lower bound of Lemma 1. We apply again Lemmas 1 and 2 to get the 
same estimate ( ^.191 ) for iP{xq) and iP^Xq + 1) at the midpoint Xq of [xi,X2] and at Xq + 1. 
By construction, L = |xo — Xg| satisfies k < L < 2(C + l)k^. Set 

, / X J ip(x), X e [xo,Xnl 
[ U, otherwise 

and 0l(x) = iPl{x)/\\iPl\\- Since Xmax e [xo,Xo], we have \\iPl\\ > |xmax| + 1- Hence 

II^lII II^lII Fmaxl + l (2.20) 

(C + 2)e-(5-2")'=/2 

for A; G iT, fc > Kq, where Kq^e, D, d) is sufficiently large. 

By the variational principle, there exists a point E' in the spectrum of Hai^i (here 
6' = ia ~ aOxinax + 6) such that 



|E'-E|< ||(iJ,,,,, -E)0i||< 

\\{H^.^e'-H^,e)(t)L\\ + ||(//a,e - ^o)0l|| + |^ - ^o| < 
C'\a - a'\L + 4(C + 2)e-(^-2e)fc/2 ^ g-fcs/4 ^ 

C"2(C + l)|a - a'|A;3 + 4(C + 3)e-»'=/^ 



-(9-2e)fc/2 , -kg/i . (2-21) 



where we applied the estimate \f{an + 9) — f{a'n + 9')\ < C'\a — a'\L with some C > and 
used that, by our choice of parameters, 2e < g/2. Let a' be sufficiently close to a, so that 
k = [4 1 In |a — a' 1 1/(7] (or /c ± 1) is in i(' and larger than Kq and | In |q; — q;'| | is sufficiently 
large depending on values of C, C, and g. Then we obtain from (|2.21|) the statement of 
Theorem 3 with c(a, 6) = 1^C\C + 1)/^^ ^ 1. ■ 



3 Proof of Theorem 1 

As we noted in the introduction, the theorem for a ^ fi is easy to prove following P, |16 . 



Take a G fi and consider the sequence of its canonical rational approximants p„/g„. Because 
of continuity in 6^, the set S'(p„/g„) consists of at most g„ disjoint intervals, say S{pn/<in) = 

n' 

UjJ\[aj,6i], g'„ < Qn- Given the continuity result of Theorem 3, the proof of Theorem 1 
differs from the corresponding analysis for the almost Mathieu operator in |T^ only in one 
detail. By [|l^, |17| we have |5'(a)| > \misxi\in^^\S{pnl<in)\- We are going to show that 
15(0;) I < liminf^^oo \S{pn/<in)\ when 7(-E) > for a.e. E. For all n sufficiently large (such 
that \pn/<in — «| < h{a, 5)), Theorem 3 says that, except perhaps for the points belonging to 
A{5), all the other points in the spectrum of H^ are necessarily inside the set (a' = Pn/^ln) 



Therefore, 



Ui=i[ai — c\{a — a') \v? \a — a\\,hi + c\{a — a') \t? \a — a'\ 



\S{a)\ < \S{pn/qn)\ +2cg„|(a-p„/gn)ln^ \oi-Pn/qn\\ +5 (3.22) 



As n — i> oo, the second addend at the right hand side tends to for any irrational a 
because of ( p.lO|) . The constant 5 > can be taken arbitrary small by Theorem 3. Hence, 
\S{a)\ = lim^^oo \S{Pn/qn)\- ■ 
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